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Abstract. - The dynamics of an impurity (or tracer particle) immersed in a dilute granular 
gas under uniform shear flow is investigated. A non-equilibrium phase transition is identified 
from an exact solution of the inelastic Boltzmann equation for a granular binary mixture in the 
tracer limit, where the impurity carries either a vanishing (disordered phase) or a finite (ordered 
phase) fraction of the total kinetic energy of the system. In the disordered phase, the granular 
temperature ratio (impurity "temperature" over that of the host fluid) is finite, while it diverges 
in the ordered phase. To correctly capture this extreme violation of energy equipartition, we show 
that the picture of an impurity enslaved to the host fluid is insufficient. 



Under external driving, an initially mixed macroscopic 
granular system may segregate. This ubiquitous effect is 
not always desirable in applications, and among the nu- 
merous counter-intuitive properties of granular materials, 
it has arguably been among the most studied in the last 
twenty years [HE]- Yet, it is not well understood. The rea- 
son is twofold: there is a large number of a priori relevant 
parameters involved in the description of the granular mix- 
ture; additionally, analytical results are scarce and difficult 
to obtain. Both shortcomings call for an improved fun- 
damental understanding of physical situations where the 
different physical effects at work (size asymmetry, forcing, 
collisional dissipation, etc.) can be deciphered. As a con- 
sequence, seeking for theoretical progress, we simplify the 
problem in two respects. First, we consider the impurity 
limit where a tracer particle is immersed in a driven dilute 
granular gas described by the inelastic Boltzmann equa- 
tion. Second, we focus on an inelastic version of the so- 
called Maxwell molecules [3], that can be seen as defining 
the kinetic theorist's Ising model. Indeed, the Boltzmann 
equation for elastic bodies is already notoriously difficult 
to deal with, and inclusion of collisional dissipation -an es- 
sential aspect of inter-grains collisions- leads to a far more 
complex description. The corresponding inelastic Maxwell 
model, that has witnessed an upsurge of interest in recent 
years [4j-[9] , is inspired by Maxwell's original insight [10] 
that scattering processes involving collision rates indepen- 
dent of the relative velocity of impacting particles, define 



mathematically tractable approaches [TT) . 

A valid question to raise from the outset is then that of 
the relevance of such a simplification. This issue should 
be addressed by comparison with more refined approaches. 
Among the few well accepted models accounting for the 
specifics of dissipative inter-grains collisions, the simplest 
is the inelastic hard sphere model [2j[12], where results 
pertaining to the inelastic Boltzmann equation can be 
compared to their Maxwell model counterpart. While the 
agreement concerning the Navier-Stokes transport coeffi- 
cients is only qualitative [T3], the situation improves sig- 
nificantly for sheared mixtures where the rheological prop- 
erties (shear and normal stresses) agree very well with (ap- 
proximate) analytical results and computer simulations for 
inelastic hard spheres |14j . This agreement, in a param- 
eter space of large dimensionality, gauges and establishes 
the reliability of the Maxwell model to capture important 
effects in sheared granular mixtures. We come back to the 
relevance of our approach before our concluding section, 
see below, and we further note that experiments with mag- 
netic grains also point to the usefulness of Maxwell models 

Our objective is to investigate the dynamics of an im- 
purity or tracer particle (hereafter referred to as species 
1) in an inelastic dilute gas (species 2), maintained in 
a non-equilibrium steady state by a uniform shear flow. 
This particular state is amenable to analytical treatment 
from the Boltzmann equation since it is characterized by 
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constant and uniform partial densities (rij for i = 1,2) 
as well as a uniform "granular temperature" T(t) [15] , 
while the only gradient affects the velocity profile, that 
is furthermore common to both species: u x = ay in a 
Cartesian frame \(x,y,z) in three dimensions]. This lin- 
ear velocity field defines the constant shear rate a that 
quantifies viscous heating (shearing work). This energy 
injection mechanism is balanced by dissipation in inelas- 
tic collisions, the latter being parameterized by the con- 
stant coefficients of normal restitution oiij for collisions 
between species i with j [12] . Although the total tem- 
perature T = {n\T\ + ri2T2)/(ni + n 2 ) {Ti is the partial 
temperature of species i) changes with time, the tempera- 
ture ratio T1/T2 is time independent in the hydrodynamic 
regime ( namely, for times much larger than the mean free 
time). Note that equipartition, of course, has no reason to 
hold in such non equilibrium conditions: T1/T2 7^ 1 [T7]. A 
plausible expectation is that upon taking the tracer limit 
c\ — >• (where c, = rii/ini + n 2 ) is the mole fraction of 
species i), the properties of the medium (or granular gas) 
are not affected by the presence of the tracer particles 
(or impurity). Let us consider the relative contribution 
of the tracer particles to the total energy of the system, 
E\/E = c\T\jT. Out of equilibrium, one expects that 
Ei/E actually scales like c\, and should therefore vanish 
when ci — > 0. In other words, the impurity should be en- 
slaved to the host medium, and should not change the 
macroscopic properties of the system. In this letter we 
present an example of a violation of the above expecta- 
tion: A non-equilibrium phase transition takes place, that 
allows to discriminate a disordered phase, where Xi/T in 
the tracer limit is finite -and hence E\jE = 0-, from an 
ordered phase where E\/E is finite -and therefore T\jT 
diverges. The corresponding regions of phase space will be 
worked out explicitly, which will reveal that both phases 
may exhibit unexpected reentrant features. The appari- 
tion of an ordered phase corresponds to an extreme vi- 
olation of equipartition, and has already been reported, 
but only in the case of an unforced system (a = 0) [T9"] . 
However, this work was performed at the "enslaved im- 
purity" level, neglecting the retroaction of impurity onto 
the host fluid. This precludes the derivation of the order 
parameter. Compared to the findings of Rcf. [19], we ob- 
tain explicitly the value of the transition order parameter 
Ei/ E in the ordered phase, and show that the scenario is 
more complex than reported so far, with the emergence of 
overlooked phases. Furthermore, we investigate the fate 
of the previous phases, including the new ones obtained, 
at finite shear (a 7^ 0). These new results are the most 
significant contribution of the present work. 

We consider a binary mixture of grains at low density, 
driven by a uniform shear flow. In the Lagrangian frame 
moving at the linear flow velocity u, the velocity distribu- 
tion functions /j(r, v, t) of each species become uniform, 
i.e., fi(r, v, t) = fi(V,t) where V = v u. In this frame, 
the set of coupled Boltzmann equations for the velocity 



distributions fi(V,t) read [TS] 
d d 

di fi ~ aVy ~iWj l = ? J <jW*' W 
x j 

where the Boltzmann collision operator Jij[fi, fj] is 

Jij [Vii/i,/,-] ffr,!'^'- J ,lrr ["tf/iCv'o/jC^) 

-/i(Vi)/.,-(V 2 )] . (2) 

Here, D, d = 2-K d / 2 /T(d/2) is the total solid angle in d di- 
mensions, and atij < 1 denotes the (constant) coefficient 
of restitution for collisions between particles of species i 
with j. Moreover, V[ = V 1 - (j,^ (l + a" 1 ) (a ■ g 12 )<7, 
V 2 = V 2 +/iy (1 + a^ 1 ) (CT-gi2)CT, where gi 2 = V1-V2, 
er is a unit vector directed along the centers of the two 
colliding spheres, and fiij = mi/{rrii + mj). The effective 
collision frequencies ujij for collisions i-j are independent 
of the relative velocities of the colliding particles but can 
depend on space and time through its dependence on den- 
sities rii and temperature T. They can be also seen as 
free parameters of the model. Here, since our problem 
involves a delicate tracer limit, we aim at the simplest 
possible approach: In the "plain vanilla" Maxwell model 
worked out here, the collision frequencies u)y are of the 
form bJij = iscj, where v~ Y is an effective (constant) mean 
free time. The form of ujij is closer to the original model 
of Maxwell molecules for ordinary gas mixtures [TU] . The 
plain vanilla Maxwell model has been previously consid- 
ered by several authors |19H20j in several problems per- 
taining to granular mixtures. Other choices of are 
possible, e.g. if the goal lies in mimicking as accurately as 
possible inelastic hard spheres [H], but at the price of a 
substantial increase in the complexity of the model that 
prevents to get exact results. 

In the uniform shear flow problem, the reduced shear 
rate a* = a/u is the relevant noncquilibrium parameter 
since it measures the distance from the unforced case, the 
much studied so-called homogeneous cooling state (HCS) 
[12] . Within our framework, a* does not depend on time, 
so that in general a* and a,ij are independent parameters. 
This decoupling allows us to analyze the combined effect of 
both control parameters on the dynamic properties of the 
granular mixture. This is one of the main advantages of 
the vanilla model used here in contrast to previous works 
[14] . More generally, a key feature of the Boltzmann equa- 
tion for Maxwell models is that the pressure tensor P can 
be exactly determined in terms of the shear rate and the 
parameters of the mixture (the mass ratio fi = mi/m 2 , 
the mole fraction c\ = +"-2) and the coefficients of 

restitution ay) [H1III]- The knowledge of the pressure 
tensor allows one to determine the non-Newtonian trans- 
port properties of the mixture. The energy ratio E±/E is 
also an important property that can be obtained from the 
pressure tensor. 

We now address the tracer limit (ci — > 0) for E\/E. 
After some tedious but straightforward algebra [35], an 
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interesting phenomenon is uncovered in this limit: while 
Ei / E vanishes for some parameter range -this is the naive 
expectation-, the possibility exists that E\jE ^ 0. Two 



relaxation rates, Aj° and X^' , play a key role in delimi 
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tating the two behaviours |22j : whenever A 
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the expression of E\ / E can be written as 
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an which is equivalent to neglecting the collisions among 
tracer particles themselves in the kinetic equation of f\. 

The change of behaviour (energy ratio that vanishes or 
not) is akin to an ordering process where the tracer is 
either enslaved to the host fluid {E\jE = 0), or carries 
a finite fraction of the total kinetic energy of the system 
(Ei/E 0). The latter situation, where the temperature 
ratio is divergent, will be subsequently referred to as the 
"ordered phase" , characterized by an extreme breakdown 
of energy equipartition. Loosely speaking, the system is 
invariant under the transformation c\ — > Sci in the disor- 
dered phase, where T\/T reaches a finite value in the limit 
ci — > 0. This invariance is broken in the ordered phase, 
since then T\/T cx o[ 1 . 

In light of the previous discussion, it is instructive to 
analyze some special cases. First, for sheared elastic 
gases (a* 7^ and ai 2 — 0122 = 1), it is easy to see 
that A 2 °^ > Aj° ("disordered" phase) if the mass ratio 
ji > \J~2 — 1 ~ 0.414 for any value of a*. However, if 
/i < \[2 — 1, X^ > A 2 ("ordered" phase) for a* larger 
than a critical value a*([i). We recover here the transi- 
tion already found |23j for elastic Maxwell mixtures under 
uniform shear flow. A second interesting situation is that 

(5) 

of the HCS (a* = but ay 7^ 1). In this case, there 
are also two different phases depending on the relative 



(4) 



where the dependence of mi/m 2 , Oy, and a* is implicitly positions of A^ 0) and X 2 V ' . In particular, the transition 
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assumed on the right-hand side. In addition, A^ 1 is de- 
fined by the expansion Ai(a*,ci) « \ f\a*) + X[ 1 \a*)c 1 + 
0{c\) where Ai is the largest real root of a sixth-degree 
polynomial equation [35]. The general expressions of D, 
A01 and Ai are too lengthy to be written down here. For 
the sake of illustration, we give their expressions for the 



point (A^ = X 2 V> ) -where it can be shown that the sys- 
tem relaxation time diverges- leads to the critical mass 
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with /ijj CS < Mhcs- The ordered phase (A^ 
exists for mass ratios smaller than /ijj CS or larger than 
/ijjog. Moreover, the first ordered phase (/j, < Hj[cs 
corresponding to light impurities) is present only when 
ai2 > + a 22 )/2, and so is absent when a.12 = 0*22 
or when ai 2 < l/y/2. It therefore appears that the occur- 
rence of this ordered phase not only requires collisional dis- 
sipation but more precisely, asymmetric dissipation. In 
the absence of shear (a* — 0), the order parameter E\jE 
can be cast in a simple form as 
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It must be noted that in general the expression of Ei/E 
does not explicitly depend on the coefficient of restitution 



The existence of the second ordered phase (/x > /%cs> 
heavy impurities) was already found by Ben-Nairn and 
Krapivsky [T!5] in their analysis on the velocity statistics 
of an impurity immersed in a uniform granular fluid. It 
must be remarked that a similar extreme breakdown of the 
energy equipartition has also been reported in the HCS 
for inelastic hard spheres [24] since, in the ordered phase, 
the ratio of the mean square velocities for the impurity 
and fluid particles T\ m 2 /T 2 toi is finite for extremely large 
mass ratios (mi/m 2 — > 00). 

We now turn to the general case where viscous heating 
and collisional dissipation are both at work. Unlike with 
the HCS where the light impurity ordering requires asym- 
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Fig. 1: Phase diagram in the reduced shear rate vs mass ratio 
plane, for different coefficients of restitution. The lines display 
the boundaries between ordered and disordered phases; the 
light and heavy impurity ordered phases are indicated. The 
vertical dot-dashed line corresponds to a cut leading to the 
inset of Fig. [2] The arrows shows Mhcs an d Mth f° r a i2 = 0.9 
and «22 = 0.6. For the two other parameters sets, {ai2, 022} = 
{0.9, 0.9} and {0.6, 0.9}, the thresholds fj,^ and /%cs ex i s t ( n °t 
shown, for clarity) while /ijjQg is not defined. Here, d — 2, but 
very similar results are obtained in three dimensions. 



metric dissipation, we find here that adding the "shear 
dimension" to the phase diagram ensures the systematic 
existence of the light impurity ordered phase. This is il- 
lustrated in Fig. [TJ see the left hand side. On the other 
hand, the heavy impurity order sets in at large enough 
mass ratio and small enough shear (lower right corner of 
Fig. [1]). As anticipated from the study of the unforced 
system, the behaviour differs depending if the inequality 
012 > v0~+ a 22)/2 holds or not. In Fig. [TJ the param- 
eter set fulfilling the above constraint is ai% = 0.9 and 
ct22 = 0.6. In the low shear rate limit, the figure indeed 
exhibits two distinct ordered phases, one at low mass ra- 
tio /x, and one in the opposite heavy intruder limit. This 
behaviour is fully consistent with the previously analyzed 
HCS scenario. Conversely, when a 12 < \/JT+ « 22 
the light impurity ordered pocket is more restricted, and 
confined in a portion of high shear, a* > a*, with quite 
a complex dependence on the coefficients of restitution 
«22 and a 12 (see fig. [TJ). In all cases, the light im- 
purity phase can only exist provided \i does not exceed 
some dissipation dependent threshold, \i < fi t h, with 
fith = v2(l + Q!i2)/(1 + 0:22) — 1- A s can be seen in 
Fig. [TJ the order/disorder threshold shear rate diverges as 



fi -> /ith- Of course, jU t h 



1 for elastic gases, as dis- 



cussed above. The phase diagram shown in Fig. [JJ also in- 
dicates that the ordered phase exhibits reentrant features. 
Upon increasing mass ratio at fixed reduced shear rate, or 
conversely increasing shear rate at fixed mass ratio, the 
following sequence may be observed for some portion of 
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Fig. 2: Order parameter of the non equilibrium transition, as 
a function of mass ratio p = mi/mj in the unforced case for 
a two-dimensional system. Here, «i2 = 0.9 and 0-22 = 0.6 The 
inset shows the same quantity as a function of reduced shear 
rate, for p = 0.03 < Mhcs- The hatched regions indicate the 
ordered phases (reentrant light impurity phase in the inset). 
The two thresholds in the unforced case are indicated by the 
vertical arrows (/%cs — 0-041 and a'hcs — 9.83). 



phase space: a first transition from ordered to disordered 
states, followed by a reverse disorder — > order transition. 
This is illustrated in Fig. [2J where the fraction of the to- 
tal energy that is transported by the impurity is plotted 
against either mass ratio, or shear rate. It is apparent 
that, for asymptotically large shear rates, the tracer con- 
tribution to the total energy can be even larger than that 
of the excess component. 

Before concluding, we come back to the reliability of our 
"plain vanilla" Maxwell approach. Although approximate, 
it turns out to provide a relevant framework for granular 
binary mixtures under uniform shear flow. This is illus- 
trated in fig. [3] for an equimolar mixture in the steady 
uniform shear flow state, namely, when viscous heating 
and collisional cooling cancel each other and consequently, 
the (reduced) shear rate a* depends on the coefficients of 
restitution ctij . Figure [3] compares the predictions of the 



plain inelastic Maxwell model to the Monte Carlo results 
obtained within the more realistic inelastic hard sphere 
model. It can be seen that the plain Maxwell model pre- 
diction captures the important trend observed within a 
more refined framework ( the discrepancies between theory 
and simulation are less than 15%, even for strong values 
of dissipation). We therefore expect that the transition 
found in this paper is not artefactual but is a robust fea- 
ture. 

In conclusion, we have analysed the tracer limit of an ex- 
act solution [14|l21j of the Boltzmann equation for a binary 
granular mixture of inelastic Maxwell gases. This solution 
applies to arbitrary values of the shear rate a and the pa- 
rameters of the mixture, namely, the mole fraction ci, the 
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Fig. 3: Temperature ratio T1/T2 versus the (common) coef- 
ficient of restitution [a = an — an = Q22) in the steady 
uniform shear flow state. The predictions of the inelastic 
Maxwell model (continuous curve, present work) are tested 
against Monte Carlo data (symbols, from [25] ) . The param- 
eters are ci = 0.5 (equimolar mixture) for a three dimensional 
system (d = 3). Two mass ratios are considered, corresponding 
to moderately heavy intruders. 



mass ratio /z, and the coefficients of restitution aij. We 
have argued that when the granular system is driven by an 
externally imposed (uniform) shear, the inelastic Maxwell 
model correctly encodes the physics of more complex mod- 
els, as evidenced by the dependence of the non-Newtonian 
transport properties on the numerous parameters of the 
mixture |14j . 

Within the framework of the mean-field Maxwell model, 
we have shown a quite unexpected result: the relative con- 
tribution of the tracer species (or impurity) to the total 
properties of the mixture does not necessarily tend to zero 
as ci — > 0. Consequently, the seemingly natural "enslaved 
impurity" picture |26] breaks down. This surprising result 
extends to granular gases some related results derived time 
ago for ordinary sheared mixtures [23j . The phenomenon 
discovered here, arising from an exact solution of the set 
of coupled Boltzmann equations, has been illustrated with 
the energy ratio E\ JE as a probe, but identical conclusions 
can be drawn for other properties of the system such as 
the (intrinsic) shear viscosity and the viscometric func- 
tions [27]. 

The corresponding extreme kinetic energy equipartition 
violation can be seen as an ordering transition, governed 
by the competition between two characteristic relaxation 
frequencies, and . We have found two different 
classes of ordered phases (see the shear rate versus mass 
ratio diagram in fig. [I}: a light impurity phase that exists 
when /i < /zjjQg for shear rates larger than a certain critical 

value and a heavy impurity phase that requires fi > (J^Qg 
and shear rates smaller than a certain threshold value. As 



fig- HI clearly shows, both ordered phases exhibit reentrant 
features. 

Usually, in the tracer limit, one assumes that the state 
of the excess component 2 (or granular gas) is not dis- 
turbed by collisions with tracer particles 1 (or impurity) 
and that the self-collisions among tracer particles can be 
also neglected [25] . The results derived in this letter show 
that while the second hypothesis is correct, the first expec- 
tation fails in the ordered phase. In this sense, the tracer 
limit must be taken with care since it presents, in the 
ordered phase, a similar complexity as the general prob- 
lem with finite concentration. In addition, although our 
present description has focused on the energy ratio for 
a binary system, we expect that the results reported here 
should have interesting consequences for polydispcrse mix- 
tures, or in terms of the velocity statistics, not addressed 
here. Finally, an important simplification allowed by our 
modelisation lies in the decoupling between shear rate and 
dissipation, while granular gases exhibit -in the steady 
state- inherent coupling between inelasticity and spatial 
gradients. This means that a fingerprint of our scenario 
should be sought for experimentally in the hydrodynamic 
transient regime, before the steady state where collisional 
dissipation and viscous heating balance each other. We 
hope that this letter stimulates the performance of such 
experiments and/or computer simulations to detect the 
transition phenomenon reported here. 
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